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Abstract 

The impact of uniaxial and hydrostatic stress on the ballistic thermal conductance (m) and the 
specific heat (C v ) of [100] and [110] Si nanowires are explored using a Modified Valence Force 
Field phonon model. An anisotropic behavior of ki and isotropic nature of C v under strain are 
predicted for the two wire orientations. Compressive (tensile) strain decreases (increases) C v . The 
C v trend with strain is controlled by the high energy phonon sub-bands. Dominant contribution of 



the low/mid (low/high) energy bands in [100] ([110]) wire and their variation under strain governs 
""O the behavior of Ki. 
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The wide application of Silicon nanowires (SiNWs) in areas ranging from Complementary- 
Metal-Oxide-Semiconductor (CMOS) transistors [U [2] to thermoelectric (TE) devices [3], 
non-volatile memories |Ij and, solar cells [3] has made SiNWs an extremely useful nano- 
electronics device concept. Device properties can be engineered [HI [7] through the introduc- 
tion of built-in stress by process engineering and by external forces like wafer bending [8] . 
Tuning thermal properties using strain can be beneficial for cooling and increasing gain of 
lasers jH |9] or improving the efficiency of TE devices [9] . 

Prior Work: Strain/stress effects on the thermal conductivity of doped bulk zinc-blende 
(ZB) semiconductors at low temperatures has been well studied [TUl [Hj. Most of the bulk 
studies in ZB semiconductors found decrease in thermal conductivity with tensile strain 
which is qualitatively attributed to (i) a decrease in the phonon mean free path [12], (ii) 
an increase in the Debye temperature [U [12], and (iii) a change in the material stiffness 
[15] . The few experimental [Hj and theoretical [HJ, [9j [131 CS] efforts that have focused on the 
strain effects on the thermal properties in nanostructures, are either limited to single crystal 
orientation (mostly [100]) or diffusive phonon transport even in small nanostructures. This 
work investigates the external strain effects on the thermal properties in ultra-scaled SiNWs 
in the coherent phonon transport regime where the wire cross-section sizes are comparable 
to the phonon wavelengths [16] (X p h = hV sn d/kBT SO oK ~ lnm at V sn d = 6.5 km/sec for a 
3nm X 3nm [100] SiNW [17], where k% and h are the Boltzmann constant and Planck's 
constant, respectively). 

Free-standing SiNWs' phonon dispersion are studied using the modified valence force 
(MVFF) model [T71 - [l9] . The MVFF model has been shown to correctly reproduce the 
strain effects like the Gruneisen parameters [181 H9] and third order elastic constants in bulk 
Si and Ge [18] . The same set of Si MVFF parameters [181 119] are used in this work to model 
the SiNW phonon dispersion under hydrostatic and uniaxial stress. 

The ballistic thermal conductance (ki) across a semiconductor wire maintained under a 
small temperature gradient AT, can be calculated from the calculated phonon dispersion, 
using the Landauer's formula [20] using Eq.Q [211 [22], 



K '( e ' T ) = J 4(z,T), where 
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where N is the number of energy bins in the entire phonon dispersion energy range. The 
terms e, M(E), T, and e are the strain percentage, the number of modes at phonon energy 
E, the temperature, and the electronic charge, respectively. The term k](E,T) is the contri- 
bution to the total Ki of SiNW for E iyrnin < E < E itmax . This energy resolved representation 
in Eq. Q allows for better understanding of the variation in Ki due to strain. 

From the calculated phonon dispersion the constant volume specific heat (C v ) at a given 
T, can be calculated (23] as, 



where m uc is the mass of the SiNW unitcell in kg. The quantity E nA is the phonon eigen 
energy associated with the branch 'n' and crystal momentum vector £ q'. 

SiNW details: Square NWs with width (W) and height (H) = 3nm with two channel 
orientations of [100] and [110] are studied (inset of Fig. [l^,b). Hydrostatic deformation 
(equal deformation along all the directions) and uniaxial deformation (along the wire axis) 
are applied to these wires varying up to ±2%. The outer surface atoms in these NWs are 
allowed to vibrate freely. Extremely small SiNWs may show significant surface and internal 
atomic reconstruction [91 [23] or phase change [9] under strain leading to larger changes in 
C v and Ki, are outside the scope of the present study. 

Ballistic thermal conductance Ki : In all the SiNWs Ki is calculated using Eq.Q at 300K. 
K\ increases (decreases) under uniaxial compression (tension) for both the wire orientations. 
Similar variations in K\ for [100] SiNWs are also observed by Li. et al(9] using non-equilibrium 
molecular dynamics (NEMD) calculations. The [100] SiNW shows larger variation in k\ 
under tensile uniaxial deformation compared to the [110] SiNW (Fig. [l^,b). Ki shows a 
weak hydrostatic stress dependence in [100] SiNW in contrast to the [110] SiNW which 
shows a decrease (increase) of 2.9% (1.98%) in K\ under 2% compressive (tensile) strain 
from the unstrained value (Fig. [lji,b). Similar strain behavior is obtained in SiNWs with 
W=H up to 5nm, which are not shown here for the sake of brevity. Thus, ultra-scaled SiNWs 
show an anisotropic variation in K\ under hydrostatic deformation. 

Specific heat (C v ): Variation in the phonon dispersion under strain also changes the 
C v of SiNWs. Compressive (tensile) strain decreases (increases) the C v in both types of 
SiNWs (Fig. [2k, b)). This variation in C v is similar to the bulk Si behavior as reported 
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in Refs.[l3j [15]. Under hydrostatic stress the C v decreases by ~ 2.7% at 2% compression 
and increases to ~ 2.6% at 2% expansion compared to the unstrained C v value for both 
the SiNW orientations. Under uniaxial deformation the variation is < 1% for ±2% strain 
for both wire orientations (Fig. [2^i,b). Hence, C v variation is isotropic under strain in 
ultra-scaled SiNWs. 

Reason for the variation in K\ : To understand the variation in Ki with strain, the energy 
dependent contributions are analyzed using Eq. Q . The entire phonon spectrum is grouped 
into three energy ranges, (i) El ('low' bands), (ii) Em ('mid' bands) and, (iii) Eh ('high' 
bands) (see Fig. |3]a). These bands show variable contributions under strain which determine 
the overall effect on Ki (Table [T]). 

Under uniaxial strain, in [100] SiNW contribution to Ki is mainly from El + Em, which 
decreases from ~88% under compression (Cm) to ~84% under tension (Tn) (Fig. [3)3), 
whereas in [110] SiNW the contribution to k\ is mainly from El + Eh, which reduces 
from ~59% (Cm) to ~54% (Tn) (Fig. gl). Under hydrostatic strain, in [100] SiNW main 
contribution to k\ is from Em, which reduces from ~56% (Cm) to ~45% (Tn) (Fig. [3jj), 
while in [110] SiNW main contribution to Ki is from El + Eh, which increases from ~52% 
(Cm) to ~58% (Fig. [3^). These details are also summarized in Tablejlj Thus, low/mid energy 
contribution for [100] SiNW and low/high energy contribution for [110] SiNW explains the 
anisotropic nature of K\. 

Reason for C v variation: Under the action of strain the contribution of the phonon bands 
in Eh range varies considerably as shown in Fig. |4j Under uniaxial stress, both [100] and 
[110] SiNWs show almost no variation in the El and Em range, and minute variation in 



TABLE I. Value of band contribution to k; under strain (Fig. |3]). Variation shown from compression 
(Cm) to tension (Tn) for uniaxial (U) and hydrostatic (H) strain for both SiNW orientations (Or). 
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the Eh range(Fig. |4^l,c). The contribution increases from Cm to Tn in the Eh range which 
governs the C v trend under uniaxial strain (Fig. [2]). Under hydrostatic stress, contribution 
from the Eh range increases from Cm to Tn, but in a much larger magnitude compared to 
uniaxial stress, for both the wire orientations (Fig. |4|D,d). This explains the larger variation 
in C v under hydrostatic stress (Fig. [2j. Thus, the higher energy bands decide the strain 
behavior of C v under strain, for ultra-scaled SiNWs. 

Conclusions and Outlook: The impact of strain on the thermal properties of ultra-scaled 
SiNWs has been provided. The ballistic thermal conductance behaves an-isotropically under 
strain, however, C v shows an isotropic behavior. The observed trends in the thermal prop- 
erties can be understood by the different types of contribution of phonon modes in different 
energy ranges. The behavior of C v in SiNWs is similar to bulk Si however, Ki variation is 
very different from bulk Si. Under the low stress condition hydrostatic strain can be ben- 
eficial in engineering C v and uniaxial stress for engineering ki for cooling lasers. However, 
high external strain (|e| > 2% ) field may be needed to engineer k\ to improve the efficiency 
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FIG. 1. Variation in k\ with hydrostatic and uniaxial stress in 3nm x 3nm SiNW with (a) [100] 
and (b) [110] orientation. The uniaxial stress is applied along the wire axis. Insets also show the 
atomic structures of both the SiNWs with the coordinate axes. 



770 

760 

750 

^ 740 
ai 

730 

^3 

u 720 
710 
700 



-o-Uniaxial 


■V- Hydrostatic 2.31% 












SiNW [100] 


-2.7% 


3nm X 3nm 


(a) 


T=300K 



760 
750 
740 
730 
720 
710 



-o-Uniaxial 

■V- Hydrostatic 2.6% 
0.55% 

■0.68% 

'SiNW [110] 
3nm X 3nm 

-2.6% T=300K 

(b) 




-1 1 

e (%) 



1 2 

e„ (%) 



FIG. 2. Variation in the C v with stress in 3nm x 3nm SiNW with (a) [100] and (b) [110] orientation. 
Two types of stress are applied in these SiNWs, (i) hydrostatic pressure and, (ii) uniaxial stress 
along the wire axis. 
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FIG. 3. (a) Phonon dispersion of 3nm x 3nm [100] SiNW with the three energy ranges considered 
for ki analysis. The percentage contribution to m from different energy ranges, for [100] SiNW 
under (b) uniaxial stress, (c) hydrostatic stress and, for [110] SiNW under (d) uniaxial stress and 
(e) hydrostatic stress. Arrows show the variation in contribution from compression (Cm) to tension 
(Tn). 
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FIG. 4. Variation in CV with phonon energy bands under various conditions. 3nm x 3nm [100] 
SiNW with (a) uniaxial stress, (b) hydrostatic pressure. 3nm x 3nm [110] SiNW with (c) uniaxial 
stress, (d) hydrostatic pressure. The higher sub-bands show larger variation in CV contribution 
compared to the lower energy sub-bands. 
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